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ABSTRACT 

Globular clusters (GCs) are emerging as a new class of 7-ray emitters, thanks 
to the data obtained from the Fermi Gamma-ray Space Telescope. By now, 
eight GCs are known to emit 7-rays at energies >100 MeV. Based on the stellar 
encounter rate of the GCs, we identify potential 7-ray emitting GCs out of all 
known GCs that have not been studied in details before. In this paper, we report 
the discovery of a number of new 7-ray GCs: Liller 1, NGC 6624, and NGC 6752, 
and evidence for 7-ray emission from M80, NGC 6139, and NGC 6541, in which 
7-rays were found within the GC tidal radius. With one of the highest metallicity 
among all GCs in the Milky Way, the 7-ray luminosity of Liller 1 is found to be 
the highest of all known 7-ray GCs. In addition, we confirm a previous report 
of significant 7-ray emitting region next to NGC 6441. We briefly discuss the 
observed offset of 7-rays from some GC cores. The increasing number of known 
7-ray GCs at distances out to ~10 kpc is important for us to understand the 
7-ray emitting mechanism and provides an alternative probe to the underlying 
millisecond pulsar populations of the GCs. 

Subject headings: Gamma rays: stars — Globular clusters: individual (Liller 1, 
M80, NGC 6139, NGC 6441, NGC 6541, NGC 6624, and NGC 6752) — Pulsars: 
general 
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Introduction 



Globular clusters (GCs) are very efficient in forming close binary systems and their dece- 
dents such as low-mass X-ray binaries (LMXBs), cataclysmic variables, and millisecond pul- 
sars (MSPs). It is commonly believed that the high stellar encounter rate o f GCs facilitates 



the f or mation of these bili ary s ystems through dyna mical interactions (e.g., IVerbunt &: Hut 
19871 ) . iPooley et al.l ( l2003l ) and iPooley fc HutI (120061 ) have shown by X-ray observations that 
the number of clo se binary s ystern s in a GC scales with the stellar encounter rate (Fc) of 
the GC. Recently, iHui et al.l (120101 ) found that a cluster with higher Fc and metallicity hosts 
more MSPs. 

Radio and X-ray observ ations have revealed about 140 MSPs in 26 globular clus- 
tersQ (IManchester et al.ll2005l ). However, the pr esence of much stron ger X-ray emitters can 
contaminate the X-ray observations of MSPs. lAbdo et al.l (l2009bl ) reported a number of 
MSPs in the Galactic field using data obtained using the Fermi Gamma-ray telescope. Be- 
cause MSPs are the only known steady 7-ray sources in GCs, 7-ray observations serve as an 
alternative channel in studying the underlying MSP populations in GCs. 

The Large Area Telescope (LAT), onboard the Fermi Gamma-ray Space Telescope, is 
a pair-production telescope designed to detect 7-rays with energies between ~20 MeV and 
>300 GeV. It operates in a survey mode in which it scans the whole sky every 3 hours. The 
68% containment radius of individual photons is 0?6 at 1 GeV, the typical energy of the 
photons collected from the GCs. The point-source sensitiv ity of LAT is ;^10~^ phcm~^ s~^ 

2009h Such sen- 
a), including 



above 100 MeV in one year of survey-mode observations ( Atwood et a. 



sitivity has enabled the discove ry of 7-rays fro m 8 GCs (Abdo et al. 



47 Tucanae ( lAbdo et al.l l2009al ) and Terzan 5 (iKong et al 



201 



2010 



a 



From the known 7-ray 



luminosity of individual MSPs in the Galactic field ( lAbdo et al.ll2009bl ). they would not be 
detected at distances of several kpc. It is therefore believed that the 7-rays from GCs do 
not come from individual MSP, but from a population of MSPs. 



The ra diation mechanism of 7- rays remains unclear, 
model 



le.g., 



In the pulsar magnetosphere 



Venter fc de Jagerll2008[ ). 7-rays up to a few GeV is radiated from the MSPs 



through curvature radiation. On the other hand, inverse Compton (IC) processes resulted 
from energetic particles up-scattering low-energy p hotons, such as starlight and infrared 



light, may give rise to 7-rays of MeV to TeV energies ( iBednarek fc Sitarekil2007l : ICheng et al. 



20101). It is worth noting that 7-ray observations at >10 GeV have not yet resulted in any 
detections (lAnderhub et al.ll2009l : lAharonian et al.ll2009[ ). 



^http: / /www. naic.edu/~pfreire/GCpsr. html 
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In this work, we searched for 7-rays from those GCs with high encounter rate. We report 
the discovery of 7-rays from the GCs Liller 1, M80, NGC 6139, NGC 6541, and NGC 6752. 



2. Cluster sample and Fermi Observations 



It has long been suggested that the stellar encoun ter rate of a GC is a measure of the 



numb er of LMXBs, the precedents of MSPs (see, e.g. IVerbunt fc HutI llQSTt iGendre et al. 



20031 ). Assuming that each MSP emits similar amount of 7-rays, the expected 7-ray lumi- 
nosity of a GC is proportional to the number o f MSPs it h o sts, an d thus to the encounter 
rate Fc, as supported by the findings made by lAbdo et al.l (j2010d ). Here we estimate the 
encounter rate by Fc oc plr^/cr where po is the central luminosity density, Vc the core radius, 
an d a the velocity disp ersion at the cluster center. We adopted th e values of a as presented 



m 



Gnedin et al.l ( 12002| ). and those of po and Vc in iHarrid (1 19961 . 2003 version). Then we 



ra nked the GCs according to the encounter rate divided by the squared distance compiled 



m 



Harris! (119961 . 2003 version), giving a measure of the relative expected 7-ray flux for each 



GC. Most of the eight known 7-ray GCs (47 Tucanae, u Centauri, M62, NGC 6388, Terzan 5, 
NGC 6440, M28, and NGC 6652) and three source candidates (NGC 6541, NGC 6752, and 
M15) are ranked high, i.e. at the top 23 in the ranking. Only NGC 6652 gave a relatively 
low expected 7-ray flux. Among the top 20, we identify the following GCs that have not 
been reported as 7-ray source or source candidate before (in descending order of expected 
7-ray flux): Liller 1, M22, NGC 2808, NGC 362, NGC 6540, NGC 1851, Terzan 6, M80, and 
NGC 6397. On the other hand , nearby 7-ray emission has been reported for NGC 6441 and 
NGC 6624 (lAbdo et al.ll2010d ). In this work, we searched for 7-rays from the GCs in the 
above list, as well as looked deeper into the cases of NGC 6441, NGC 6541, NGC 6624, and 
NGC 6752. NGC 6139 (ranked 30th) was also studied. 

The 7-ray data used in this work were obtained between 2008 August 4 and 2010 
August 21. We used the Fermi Science Tools v9rl5p2 package to reduce and analyze the 
data provided by the Fermi Science Support CenteiB. Only those data that passed the 
most stringent photon selection criteria (i.e. the "diffuse" class) were used. To reduce the 
contamination from Earth albedo 7-rays, we excluded events with zenith angles greater than 
105°. The instrument response functions (IRFs) "P6_V3_DIFFUSE" were used. 



^http:/ /fermi. gsfc.nasa.gov/ssc/data/analysis/scitools/ 
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3. data analysis and results 



Photons with energies between -Emin and -Emax that come from a circular region-of- 
interest (ROI) around the GCs were included in the analysis. The ROl, as well as -Emin 
and -Emax; for those GC coincident with 7-ray emission are given in Table [TJ We sub- 
tracted the backgro und contribution from all sources in the first Fermi/LAT catalog (IFGL; 
Abdo et al.l l2010bl ) within the circular region of 15° radius around the GC position, as 
well as diffuse emission by including Galactic diffuse model (glLiem_v02.fit) and isotropic 
background (isotropic_iem_v02.txt). We assumed single power laws for all first Fermi/LAT 
catalog sources considere d, except for 7-ray pulsars of which the spectra follow power laws 
with exponential cut-off (lAbdo et al.ll2010al Jd). 



We created the "Test-statistic (TS: lMattox et al.lll996l ) maps" in the neighboring region 
of the clusters using the tool gttsmap. These TS maps are created by moving a putative 
point source through a grid of locations on the sky and maximizing — log(likelihood) at each 
grid point, while stronger and presumably well-identified sources are included in each fit. 
Weaker sources would be identified in these maps. Such TS maps around Liller 1, M80, 
NGC 6139, NGC 6441, NGC 6541, NGC 6624, and NGC 6752 are shown in Figs. □ to S 
centered on the best-fit centroid of the 7-ray emission (to be determined in the following). 
Photons with energy between i?min and -Emax as stated in Table [T] were used in generating 
the maps, except for Liller 1 whose map was generated using photons with energy between 
2-30 GeV to avoid contamination from the strong Galactic background due to its Galactic 
position of / = 354?84 and 6 = -0? 16. 

We found evidence for 7-ray emission within the tidal radius of Liller 1, M80, NGC 6139, 
NGC 6541, NGC 6624, and NGC 675 2. In addition, we fo und significant 7-ray emission next 
to NGC 6441, previously reported in lAbdo et al.l ( l2010d ). In all these cases, a 7-ray source 
candidate in the neighborhood of the corresponding GCs is detected above the modeled 
background emission. We did not find significant 7-ray emission within the tidal radius 
of other candidates: M22, NGC 2808, NGC 362, NGC 6540, NGC 1851, Terzan 6, and 
NGC 6397. 

The TS values and the corresponding significances (= y/TS) of the 7-ray emission 
coincident with the GCs are shown in Table [2j The number of gamma-ray photons associated 
with the corresponding 7-ray source in the best-fit model provided by the likelihood tool gtlike 
are also listed. Since the significance values for a number of positions were calculated, this 
significance has to be corrected according to a "trial factor" to give a "post-trial" significance. 
This trial factor accounts for the increased probability of finding a fake signal with a number 
of search positions. We estimate the trial factor as A^gc x A'bin, where A'gc is the number 
of GCs searched, and Abin the number of bins within the tidal radius of each GC. We have 
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searched for 7-rays from 30 GCs(including th ose with the highes t encounter rate normahzed 



by its distance and those first mentioned in lAbdo et al.ll2010q ). therefore A^^gc = 30. For 



Abin, the number of bins is obtained by dividing the area of the tidal radius of an GC (we 
take the tidal radius of Liller 1, 13', as the averaged tidal radius) by the area of each bin 
(0?1 X 0?1 = 0.01 squared degrees). This gives Nun ~ 15. After accounting for the trial factor, 
the "post-trial" detection significance values are obtained. Since the TS value for adjacent 
grid points are correlated in the likelihood analysis procedure, the post-trial significance can 
be considered as conservative estimates. 

We then used gtfindsrc to determine the positions of the 7-ray source candidates and 
gtlike (unbinned likelihood analysis) to obtain their spectra and their TS values. The inte- 
grated photon fiuxes and energy fiuxes are given for the energy range 100 MeV to 100 GeV, 
using simple extrapolation from the respective models. For the three cases of most significant 
detection (i.e., Liller 1, NGC 6624, and the 7-ray source next to NGC 6441), photons be- 
tween 300 MeV and 100 GeV are divided into 6 logarithmically-equally spaced energy bins. 
The fiux in each bin was obtained by gtlike using photons in the corresponding bin, while 
both the normalization and photon indices of all 7-ray sources and the diffuse background 
components, namely Galactic diffuse and isotropic components, were set to be free. Only for 
cases where the photon statistics is not enough for a satisfactory fit, the photon index was 
fixed at the best-fit value shown in Table [2l 

We also performed a long-term temporal analysis of the 7 7-ray sources, in which the 2- 
year data were binned into several periods depending on the photon statistics. No significant 
7-ray variability nor fiaring period was found, indicating that the sources are stable in 
radiating 7-rays. 

In the following section, we discuss the results of each of the GCs in more details. 



4. Gamma-ray properties from individual clusters 
4.1. Gamma-ray emission within cluster's tidal radius 

4.1.1. Liller 1 

The encounter rate of Liller 1 is the highest among the target list mentioned in Section |5J 
As seen in Fig. [H 7-ray emission from the direction of Liller 1 was found. Fitting the 
photon spectrum (dN/dE) with the single power law, the TS value of this source is 107, 
corresponding to a detection significance of 10.3cr (i.e. 9.4cr post-trial). The best-fit position 
is R.A. = 263?20, Dec. = -33?39 ±0?04(stat) (J2000; 1-a uncertainty). For localization of 
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a typical 7-ray s ource, the systema tic position error is estimated to be about 40% of the 



statistical error jAbdo et alJboQcl lN. The best-fit position is T.h from the core position of 



the GC Liller 1, which is within its tidal radius (12^57). Fig. depicts the spectrum (in 
E'^dN/dE) of the 7-ray emission between 300 MeV and 100 GeV. The photon spectrum of 
Liller 1 can be well fit by a single power law with an index of 2.2 ± 0.1, and the integral 



photon flux between 100 MeV and 100 GeV is Fq.i- 



100 GeV 



(6.9 ± 2.3) X lO-^cm-^s" 



in this model. On the other hand, the parameters provided by the power law with an 
exponential cut-off (PLE) model with all three parameters being free are not well constrained; 
we then fixed the cut-off energy in this model at 2 GeV and 10 GeV, and obtained smaller 
values of TS (i.e. 90 and 100, respectively). We therefore do not consider this model to be 
significant. We found that the rightmost spectral point shown in Fig. |5l that corresponds to 
40-100 GeV emission, were detecte d with a TS value of 16, corresponding to a significance 
of 4o". At a distanc e of 9.6 kpc (HarrisI 



(8.3±1.8) kpc given in lOrtolani et al. 



1996 



2003 version; which is within the range 

(|2Q03)), the 7-ray luminosity of Liller 1 is Lo.i-iooGeV = 
(5.9 ±2.0) X lO^^ergs"^, which is the largest among all detected 7-ray emitting GCs by now. 



4.1.2. M80 

Among the GCs in the target list, M80 was found to be a possibly 7-ray emitting GC. 
This GC is located at the Galactic coordinate of (/,6) = (352?6,19?4), i.e. well outside the 
Galactic plane. Therefore we extended the ROI radius to 10° and included photons with 
energies 200 MeV to 100 GeV in the analysis. Inside the tidal radius of M80 we found 
evidence for 7-ray emission that peaks at a TS value of 27, corresponding to a significance 
of 5.2(7 (i.e. 3.9(T post-trial), as shown in Fig. [T] (region ^4). We therefore regard M80 as a 
possible 7-ray emitter. In addition, it is possible that the possible 7-ray emission region may 
extend well outside the GC, i.e. into region B. Formally we cannot separate the whole 7- 
ray emission into two emitting regions. Assuming that the 7-ray emitting region coincident 
with the GC (i.e. region A) is a separate source, its best-fit position is R.A. = 244?23, 
Dec. = -23?02 ±0?06(stat) (J2000). Its distance of 3^3 from the GC core position makes 
the two positions consistent within errors. The photon spectrum of this 7-ray source can 
be well fit by a single power law with an index of 2.0 ± 0.2. The integral 7-ray fiux is 
-^o.i--iooGeV = (5.5 ± 3.6) X lO^^cm^^s^^. The low statistics of photons does not allow 
for accurate determination of the parameter values in the PLE model; we therefore do not 
consider this model. At a distance of 10.3;^q;7 kpc ( jBrocato et al.l Il998l ). i^o.i-iooGeV = 
8.41^1 X lO^^ergs-i. 



■^We assume similar systematic errors for other GCs as well. 
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The best-fit position of tlie 7-ray emitting region B (TS=30) is R.A. = 244?42, Dec. = 
— 23?50 ±0?07(stat) (J2000). Tfie plioton spectrum of tliis 7-ray source can be fit by a single 
power law with an index of 2.2 ± 0.2. The integral 7-ray flux is -Fo.i-iooGcV = (1-1 =t 0.6) x 
10~^cm"2s-i. 



4.1.3. NGC 6139 

The encounter rate of NGC 6139 is the 30th highest in our target list. We found evidence 
for 7-ray emission at the position of NGC 6139. The TS value of this source candidate is 
32, corresponding to a detection significance of 5.6o" (i.e. 4.5(T post-trial), obtained in the 
power law model. We therefore regard this detection as a marginal detection. The TS map 
around NGC 6139 is shown in Fig.d The best-fit position is R.A. = 246?83, Dec. = -38?90 
±0?07(stat) (J2000), which is 5!l from the core position of NGC 6139, but within its tidal 
radius (8f52). If the 7-ray source is confirmed, the 7-ray spectrum of NGC 6139 can be well 
fit by a single power law with an index of 2.1 ± 0.2. The integrated photon flux between 
100 MeV and 100 GeV is Fo.i-iooGeV = (9.9 ± 5.4) x lO^^cm-^g-^. The low statistics 
of photons does not allow for accurate determination of the parameter values in the PLE 



model; we therefore do not consider this model. At a distance of 10.1 kpc ( lHarrislll996l . 2003 



version), the 7-ray luminosity of NGC 6139 is -^vo.i_ioogcV = (1-1 ± 0.6) x lO^^ergs ^ 



4.1.4. NGC 6541 



Inside the tidal radius of NGC 6541, we found a possible 7-ray source that peaks at 
TS=19, corresponding to a significance of 4.4(T (i.e. 2.8o" post-tria l). Therefore , we co nsider 
this source to be possibly detected, consistent with the finding of lAbdo et al.l (l2010d ). The 
TS map is shown in Fig.d The best-fit position is R.A. = 272?04, Dec. = -43?85 ±0?ll(stat) 
(J2000), 8f9 from the core position of the GC NGC 6541, w e ll with in its tidal radius (29^60). 
This position is consistent with those derived in lAbdo et al.N2010bl Hl Therefore, we identify 
this source candidate to be IFGL J1807.6— 4341. A power-law fit of the photon spectrum gave 
an index of 2. 2 ±0.2 and an integrated photon flux of -Fo.i-iooGoV 
This fl ux is consistent with the upper limit of 1.1 x 10 



(9.5±4.4) X 10 -^cm-^ s"^ 



"*cm ^ s 



as given m 



AbdoetaL 



(l2010d ). noting that their flux was calculated using the PLE model. Moreover, our derived 



integrat ed energy flux of -En 1-1 nnGeV = (6.5±3.0)xl0 ^^ergcm ^s ^ are consistent with that 



given m lAbdo et al.l ( l2010bl ). The PLE model was not used in the likeli hood fit because of th e 
putative nature of the 7-ray emission. At a distance of (6.9 ± 0.7) kpc (ILee fc Carneyll2006l ). 
the 7-ray luminosity of the 7-ray source candidate is -^vo.i-iooGcV = 3.711? X 1034 ergs . 
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4.1.5. NGC 6624 



The 7-ray source IFGL J1823.4-3009 lies close to NGC 6624 flAbdo et all 120101 
Using a larger data set, we confirmed this 7-ray source and obtained a TS value of 121, 
corresponding to a detection significance of 11. Oct (i.e. 10. la post-trial). The best-fit position 
was found to be R.A. = 275?86, Dec. = -30?16 ±0?05(stat) (J2000), 12;2 from the core 
p osition of NGC 662 4 (shown in Fig. [3]). This position is consistent with that reported 



m 



Abdo et al.l (l2010d ) and it lies within the tidal radius of the cluster (20f55). We obtained 



a single power-law fit of the photon spectrum with an index 2.2±0.1, in which the integrated 
7-ray flux is -fo.i-iooGeV = (2.6 ± 0.6) X 10~^cm~^s~^. The above power-law index and our 
derived integrated energy flux o f .En.i-in nGeV = (2.0 ±0.5) x 10~^^ergcm~^ s~^ are consistent 



with that given in I Abdo et al.l (l2010bl ). On the other hand, the parameters provided by 
the PLE model (with all three parameters free) are not well constrained; we then fixed the 
cut-off energy in this model at 600 MeV, 1 GeV, 1.5 GeV, 2 GeV, 3 GeV and 5 GeV, and 
obtained TS values of 111, 120, 121, 120, 118, and 116, respectively). Fixing the cutoff 
energy at 1.5 GeV, the power law index was found to be 0.9 ± 0.2, and the integrated 7-ray 



flux was found to be -Fo.i-iooGeV = I-O.q^s x 10^ cm^ s^ . Since this model gives the same 
TS(= 121) as the power law, it is a model as good as the power-law model. We found that 
the two rightmost spectral points shown in Fig. [6l that together correspond to 14-100 GeV 
emission, were detected with a TS value of ~10 each, corresponding to a significance of 3a 



each. At a distance of 7.9 kpc (jHarrislll996l . 2003 version), the 7-ray luminosity of NGC 6624 
(using the power-law model) is Lo.i-iooGeV = (1-5 ± 0.3) x lO'^^ergs"^ 



4.1.6. NGC 6752 



We found significant 7-ray emission from NGC 6752. This source has a TS value of 49, 
corresponding to a significance of 7.0cr (i.e. 6.0a post-trial), in the power law model. As 
shown in Fig. |3l the source is compact and located close to the core of NGC 6752. The 
best-fit position is R.A. = 287?57, Dec. = -59?96 ±0?09(stat) (J2000), 4^7 from the core 
position of NGC 6752, and is consistent with being emitted from within the half-mas s radius , 
as shown in Fig. [31 This position is consistent with that derived in I Abdo et al.l (l2010d ). 
The photon spectrum of the source can be well fit by a single power law with an index 
of 2.0 ± 0.2. The integrated photon fiux between 100 MeV and 100 GeV is -Fo.i-iooGcV = 
(6.0 ± 2.8) X 10~^cm~^ s ~^. Th is fiux is consistent with the upper limit of 7 x 10~^cm~^ s~^ 
as given in lAbdo et al.l (l2010cl ). noting that their fiux was calculated using the PLE model. 



The parameters provided by the PLE model (with all three parameters free) are not well 
constrained; we then fixed the cut-off energy at 600 MeV, 1 GeV, 3 GeV and 10 GeV, and 
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obtained smaller values of TS (i.e. 23, 27, 34, and 42, respectively). We therefore d o not 
consider this model to be significant. At a distance of (4.4 ± 0.1) kpc (jHarrislll996l . 2003 
version), the 7-ray luminosity of NGC 6752 is i^o.i-iooGcV = (1-4 ± 0.7) x 10^''ergs~^ 



4.2. Gamma-ray emission outside cluster's tidal radius: the case of NGC 6441 



A 7-ray source next to NGC 6441 was first reported by lAbdo et al.l (l2010d ). Using a 
larger data set, we also found a source (a TS value of 101 is obtained with the PLE model, 
corresponding to a detection significance of 10. Oo", i.e. 9. la post-trial) with the best-fit 
position of R.A. = 267?63, Dec. = -36?89 ±0?09(stat) (J2000), 10^6 from the position of 



NGC 6 441 (shown in Fig. H]). This position is consistent with that reported in lAbdo et al. 



(l2010d ) and it likely lies outside the tidal radius of the cluster (8f00) but it is still possible 
that some emission is from within the tidal radius given the current statistical and systematic 
errors. We obtained a single power-law fit with an index 2.4±0.1, in which the integrated 7- 
ray flux is -fo.i-iooGeV = (4.1±1.0) x lO^^cm"^ s~^. This fit gives a TS value of 73. The photon 
spectrum can better be fit by a PLE where the photon index, F, is 0.4 ± 0.9 and the cutoff 
energy is (1.0±0.5) GeV, giving TS = 100. To further constrain the spectral parameters, we 
proceeded to fix F at 0.2 to 2.0 (with steps of 0.2) while letting the normalization and free, 
and found that TS=101 for F^ =0.2, 0.4, and 0.6 and the TS value decreases with increasing 
F. If one fixes F^ at 0.4, = (1.0 ±0.1) GeV is obtained. In this case, the integrated 7-ray 
fiux is -Fo.i-iooGeV = (1-0 ± 0.2) X 10~*cm~^s~^. The significance of the PLE model over 
the single power law model can be estimated by TSple — TSpi, = 28, which corresponds 
to 5.3(7. Therefore, we consider the PLE model to be a significantly better model than the 
power law model. The PLE model spectrum, as w ell as the data binn ed in different energies, 
is depicted in Fig. [71 At a distance of 11.7 kpc ( jValenti et al.ll2004l ). the 7-ray luminosity 
of NGC 6441 (using the power law model) is Lo.i_iooGeV = (3.7 ± 0.9) x lO^^ergs"^, or 
^0.1-100 Gcv = (1.3 ± 0.2) X lO^^ergs-i (using the PLE model). 



discussion 



In this work, we report the discovery of three new 7-ray emitting GCs: Liller 1, 
NGC 6624, and NGC 6752, and we found evidence for 7-ray emission from M80, NGC 6139, 
an d NGC 6541 . We a lso confirmed the 7-ray detection next to NGC 6441 first mentioned 
by lAbdo et al.l (l2010d ). We searched through various X-ray and radio catalogs for plausible 
7-ray sources in the neighborhood of the GCs. We found two pulsars next to Liller 1, as 
shown in Fig. [T] (left panel). The spin-down luminosity of PSR J1733-3322 and PSR J1734- 
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3333, is 8.3 x lO^^ergs"^ and 5.6 x lO^'^ergs'S respectively f lMorris et al.l 120021 ) . I.e. more 
than an order of magnitude lower than the observed 7-ray luminosity of lO^^ergs"^. To- 
gether with their large offsets from the 7-ray emitting region, the association of the observed 
7-rays with either pulsar is unlikely. In all cases, the GCs seem to be the only plausible 
counterparts of the observed 7-rays (we refer to region A in the case of M80). Without clear 
counterparts in other wavelengths, the nature of the 7-ray emitting region B shown in Fig. [1] 
(right panel) remains unclear, and its relation with M80 cannot be confirmed or ruled out 
at this stage. 

Several MSPs are known in NGC 6441, NGC 6624, and NGC 6752 (all located close 



to its core) but none has yet been uncovered in other cases (e.g., iHui et al.l l2009l ) . We 
speculate that there is a population of MSPs - the only known steady 7-ray emitters in 
GCs - in the 7-ray emitting GCs that we have reported here. For example, the brighter 
X-ray population of M80 is in many ways sir nilar to 47 Tucanae , the first 7-ray GC known, 
while the fainter X-ray sources may differ ( iHeinke et al.l |2003| ) . The discovery of 7-rays 
from M80, if confirmed, would suggest that it also hosts a population of MSPs. Under the 
assumption that the observed 7-ray luminosity depends solely on the number of MSPs (e.g. 
in the magnetosphe ric model), the predicted number of MSPs in the cluster using Eq.(l) of 
Abdo et al.l (l2010cl ) is listed in Table |2l Future multiwavelength observations of any MSPs 
should help in understanding the 7-ray production mechanisms and the underlying MSP 
populations. 

The 7-ray spectra are in general consistent with single power law model, and shows no 
convincing evidence for cut-off at high energies. It is mainly due to limited photon statistics. 
However, the cases of the more significant cases of Liller 1 and NGC 6624, in which 7-ray are 
found within the tidal radius, suggests that >10GeV emission exists at least for these two 
cases. If this is true, this high energy emission conform with the IC model, rather than the 
m agnetospheric model. The >10G eV emission r aay co rrespond to the IC peak as predicted 



m 



Bednarek fc Sitarekl fl2007h and ICheng et al.l f l2010l ) 



If the 7-ray emission is associated with the GCs, one striking feature is the displacement 
of the 7-ray emission from the GC cores. We briefly discuss this feature in the context of 
two radiation scenarios commonly discussed in the literature. 

In the pulsar magnetosphere model, one expects the emission to coincide spatially with 
the MSP population. While MSPs should concentrate in the GC core in general, the finite 
number of MSPs in a GC (i.e., ~10, see Tabled except for Liller 1) may invoke some offset 
of the MSPs - and thus the 7-rays they emit - from the core. Furthermore, in addition to 
the MSP population formed in the core, some MSPs are formed near the tidal radius of the 
GCs. The superposition of 7-rays from these two populations may not coincide with the 
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GC core, but is at a certain distance from the core within the tidal radius. However, this 
scenario cannot accommodate the case of the 7-ray source next to NGC 6441. 

In the inverse Compton models, the 7-ray emitting region depends not only on the 
distribution of the underlying MSPs, but also on the cooling timescale and the diffusion 
timescale of the accelerated particles. When the diffusion timescale is much smaller than the 
cooling timescale, the 7— rays may come from the outskirts of the GCs. 

In the cases of Liller 1 and NGC 6624, the gamma-ray spectrum (cf. Fig. \5\ and Fig. E]) 
can extend up to energies >40 GeV, whic h cannot be explain ed in terms of curvature ra- 
diation inside the hght cylinder. Recently, ICheng et al.l ( I2OIOI ) have suggested that inverse 
Compton scattering between background relativistic electrons/positrons in pulsar wind and 
the soft photons from the galactic disk can produce the observed gamma-rays from globular 
clusters. They predict that if the inverse Compton upscattered Galactic infrared photons 
are responsible for GeV photons, then the Compton upscattered galactic optical photons 
can produce gamma-rays up to 100 GeV. They also predict that the gamma-rays are diffuse 
and are emitted from region much beyond the core of the globular cluster. If this is true the 
center of gamma-ray emission regions is affected by three factors: the proper motion of the 
globular cluster, the asymmetric diffusion coefficient and nearby optical/IR external sources 
like stars and nearby regions containing an enhanced amount of dust. 

The discovery of several more 7-ray emitting GCs reported in this work has allowed 
population studies to be carried out in order to better u nderstand the relati onship between 
the 7-ray properties of GCs and other parameters of GCs. lAbdo et al.l (j2010d ) found that the 
stellar encounter rate scales with the 7-ray luminosity, based on a sample of 8 7-ray emitting 
GCs. While the 7-ray luminosity should correlate with the number of MSPs in the clusters in 
both models, it is expected that the low-energy photon energy density is also a good estimator 
of 7-rays in the IC models. Recently, using a total of 15 GCs that are 7-ray sour ces or source 



candid a tes, includ i ng th ose reported in this work and those previously known (lAbdo et al. 



2010d ). iHui et al.l (120 111 ) have shown that correlations also exist between 7-ray luminosity 
and metallicity, as well as 7-ray luminosity and optical/infrared photon energy density. 

To conclude, we report the detection of new 7-ray emitting GCs, based on a ranking 
list that replies on the encounter rate of the GCs. The 7-ray luminosity of Liller 1 is the 
highest of all known 7-ray GCs. Its non-detection before this work may be due to its position 
that lies on the Galactic plane, where Galactic diffuse emission is very strong. Liller 1 has 
very high encounter rate and metallicity, and may indeed host a large population of MSPs 
that have not been uncovered. Some 7-ray emitting GCs presented here are at distances of 
~10 kpc away from us, whereas lAbdo et al.l (l2010bl ) present the upper limits of known GCs 
out to a distance of 6 kpc only. Therefore, we have demonstrated that the study here has 



- 12 - 



gone a further step towards finding new sources even at larger distances. This then helps to 
expand the sample of known 7-ray GCs for further studies. 
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Tabic 1: Region-of- interest (ROI), minimum photon energy (E'min), and maximum photon energy (E'max) in the analysis 
of the GCs 



GC name 


ROI center 


ROI radius 


E ■ 


E 


R.A. (J2000) Dec. (J2000) 










n n 


n 


(GeV) 


(GeV) 


Liller 1 


263.20 -33.39 


5 


0.3 


100 


M80 


244.23 -23.02 


10 


0.2 


100 


NGC 6139 


246.83 -38.90 


6 


0.3 


100 


NGC 6441 


267.63 -36.89 


5 


0.3 


100 


NGC 6541 


272.04 -43.85 


6 


0.3 


100 


NGC 6624 


275.86 -30.16 


6 


0.3 


100 


NGC 6752 


287.57 -59.96 


10 


0.2 


100 



Table 2: Summary of 7-ray results of the GCs 



GC name 


Spectral 


TSSij 


j;nificance 




Photon index*^ 


Cutoff energ 


;yPhoton flux'^ 


Energy flux'^ 


Distance*^ 


Luminosity^ A^msp^ 




model 




(^) 




(r) 


(^c, GeV) 






(kpc) 






Liller 1 


PL 


107 


10.3 


966 


2.2 ±0.1 




69 ±23 


53 ±18 


9.6 


59 ±20 


41Q+480 
^-'-'-'-210 


M80 


PL 


27 


5.2 


125 


2.0 ±0.2 




5.5 ±3.6 


6.6 ±4.2 


10.3^7 




c-o + 123 
Oo_45 


NGC 6139 


PL 


32 


5.6 


154 


2.1 ±0.2 




9.9 ±5.4 


8.8 ±4.8 


10.1 


11±6 


7c: + 114 
' ^-50 


NGC 6441 


PL 


73 


8.5 


451 


2.4 ±0.1 




41 ± 10 


23 ±6 


11.7 


37 ±9 


h 




PLE 


101 


10.0 


380 


0.4^ 


1.0 ±0.1 


10 ±2 


8.1 ± 1.0 


11.7 


13 ±2 




NGC 6541 


PL 


19 


4.4 


109 


2.2 ±0.2 




9.5 ±4.4 


6.5 ±3.0 


6.9 ±0.7 


7+2.9 




NGC 6624 


PL 


121 


11.0 


362 


2.2 ±0.1 




26 ±6 


20 ±5 


7.9 


15 ±3 






PLE 


121 


11.0 


326 


0.9 ±0.2 


1.5 


10 ±3 


7.4 ±1.0 


7.9 


5.5 ±0.7 




NGC 6752 


PL 


49 


7.0 


130 


2.0 ±0.2 




6.0 ±2.8 


6.0 ±2.8 


4.4 ±0.1 


1.4±t).7 
1-^ 





01 



°' Modeled photon number associated with the corresponding 7-ray source in the best hkehhood fit 
All the quoted errors are statistical and Icr for one parameter of interest. 

^ Integrated 0.1-100 GeV photon flux in unit of lO^^cni^^s^^ 
Integrated 0.1-100 GeV energy flux in unit of lO^^^erg cni^^s^^ 



H 



arrijj |l99d . 2003 version), except for M80 |Brocato et al.lllQQSl ). NGC 6541 |Lee fc Carnevlbood ) . and 



Distance adopt ed from 
NGC 6441 (Val enti et aLll2004[ ). 
^ 0.1-100 GeV luminosity in unit of lO'^'^erg s~^ 

^ Predicted number of MSPs in the cluster using Eq.(l) of lAbdo et al.l pOlOd) . under the assumption that the observed 7-ray luminosity 
depends solely on the number of MSPs 
^ In this case, it is unlikely that the number of MSPs is the sole factor in determining the observed 7-ray luminosity due to the offset 
' Model parameters without quoted errors are fixed at the value given. 
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Fig. 1. — The Test-statistic (TS) maps of regions of 2°x2° centered on the best-fit centroids 
(labeled by crosses) of the 7-ray emission from Liller 1 [left) and M80 (right). The size of 
each cross indicates the 1-a statistical error in the determination of the centroid position. 
The color scale shows the TS value of every bin of an area 0?1 x 0?1. The li ght blue circles 
represent the tidal radius of the respective GCs compiled in iHarrisI (119961 2003 version). 
Left: PSR J1733-3322 and PSR J1734-3333, are marked by two white diamonds. The 68% 
error positions of the LAT first catalog sources IFGL J1732.3-3243c and IFGL J1730.0-3408c 
(treated as background) are sh own as eclipses. Right: The position of the 7-ray MSP PSR 
J1614-2230 JAbdo et al.ll2009bh is shown as a cross. 
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Fig. 2. — The TS maps of regions of 2°x2° centered on the best-fit centroids (labeled by 
crosses) of the 7-ray emission from NGC 6139 (left) and NGC 6541 ( right). The circles 
represent the tidal radius of NGC 6139 (left) and NGC 6541 (nght) jHarrisI Il996l . 2003 
version). The meanings of the color scale and the size of the crosses are the same as in 
Figure [H 
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Fig. 3. — The TS maps of regions of 2°x2° centered on the best-fit centroids (labeled by 
crosses) of the 7-ray emission from NGC 6624 (left) and NGC 6752 {right). The TS map of 
a region of 2°x2° centered on the best-fit centroid (labeled by a cross) of the 7-ray ernission 
next to NGC 6624. Left panel: The circle represents the tidal radius of NGC 6624 (IHarris 
19961 . 2003 version). The white diamond shows the position of both PSR B1820-30A and 
PSR B1820-30B. Right panel: The outer an d inner circl e represents the tidal radius and 
half-mass radius of NGC 6752, respectively (jHarrislll996l . 2003 version). The meanings of 
the color scale and the size of the crosses are the same as in Figure [H 
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Fig. 4. — The TS map of a region of 2°x2° centered on the best-fit centroid (labeled by 
a cross) of the 7-ray ern ission next to NGC 6441. The circle represents the tidal radius of 
NGC 6441 (jHarrislll996l . 2003 version). The meanings of the color scale and the size of the 
cross are the same as in Figure [H 
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Fig. 5. — The spectrum {E'^ x dN/dE) of the 7-ray emission found within the tidal radius of 
Liller 1. Photons between 300 MeV and 100 GeV are divided into 6 logarithmically-equally 
spaced energy bins and The derived flux for each bin is plotted at the corresponding mean 
energy. The solid line indicates the best-fit single power law model. The rightmost point 
has a TS value of 16, corresponding to a significance level of four. 
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Fig. 6. — The spectrum (i?^ x dN/dE) of the 7-ray emission found within the tidal radius 
of NGC 6624. Photons between 300 MeV and 100 GeV are divided into 6 logarithmically- 
equally spaced energy bins and The derived flux for each bin is plotted at the corresponding 
mean energy. The solid line and the dashed line indicate the best-fit single power law model 
and the PLE model with Ec — l.b GeV, respectively. Each of the rightmost two points has 
a TS value of ~10, corresponding to a significance level of about three each. 
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Fig. 7.— The spectrum {E^ x dN/dE) of the 7-ray emission outside the tidal radius of 
NGC 6441. The sohd hne indicates the best-fit PLE model 



